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ABSTRACT 

Recent X-ray observations of hot gas in the galaxy cluster MS 0735.6+7421 reveal huge radio-bright, 
quasi-bipolar X-ray cavities having a total energy ~ 10 62 ergs, the most energetic AGN outburst 
currently known. We investigate the evolution of this outburst with two-dimensional axisymmetric 
gasdynamical calculations in which the cavities are inflated by relativistic cosmic rays. Many key 
observational features of the cavities and associated shocks are successfully reproduced. The radial 
elongation of the cavities indicates that cosmic rays were injected into the cluster gas by a (jet) 
source moving out from the central AGN. AGN jets of this magnitude must be almost perfectly 
identically bipolar. The relativistic momentum of a single jet would cause a central AGN black hole 
of mass 10 9 M Q to recoil at ~ 6000 km s _1 , exceeding kick velocities during black hole mergers, and 
be ejected from the cluster-center galaxy. Observed deviations from bipolar symmetry in the radio 
cavities can be caused by subsonic flows in the ambient cluster gas, but reflection shocks between 
symmetric cavities are likely to be visible in deep X-ray images. When the cavity inflation is complete, 
4PV underestimates the total energy received by the cluster gas. Deviations of the cluster gas from 
hydrostatic equilibrium are most pronounced during the early cavity evolution when the integrated 
cluster mass found from the observed gas pressure gradient can have systematic errors near the 
cavities of ~10-30%. The creation of the cavity with cosmic rays generates a long-lasting global 
cluster expansion that reduces the total gas thermal energy below that received from the cavity shock 
- even this most energetic AGN event has a net cooling effect on cluster gas. One Gyr after this single 
outburst, a gas mass of ~ 6 x 10 11 Mq is transported out beyond a cluster radius of 500 kpc. Such 
post-cavity outflows can naturally produce the discrepancy observed between the cluster gas mass 
fraction and the universal baryon fraction inferred from WMAP observations. 

Subject headings: cosmic rays- cooling flows - galaxies: clusters: general - galaxies: active - X-rays: 
galaxies: clusters 



1. INTRODUCTION 

The hot gas in galaxy clusters emits prolifically in X- 
rays, and has been extensively studied by X-ray tele- 
scopes Chandra and XMM-Newton, revealing complex 
structures such as X-ray cavities, thermal filaments, and 
weak shocks in the intracluster mediu m (ICM) as in 
the Perseus cluster ()Fabian et al.l 1201 )6i. These inter- 
esting structures are presumed to be related to the ac- 
tive galactic nuclei (AGNs) located at the cluster cen- 
ters. Although the ICM in many clusters losses en- 
ergy by radiation in a timescale much shorter than the 
cluster age, high-resolution X-ray spectroscopy indicates 
that the gas does not cool to low temperatu res (e.g., 
iPeterson et"alll200lL 120031: iTamura et alll2001t for a re- 
view see IPeterson fc Fabian! I2006D . The deficit of cool- 
phase gas and young stars in cluster cores, often known 
as the 'cooling flow' problem, suggests that the ICM is 
prevented from the cooling catastrophe by one or more 
heating sources. Currently, heating of the ICM by central 
AGNs, as indicated by X-ray cavities and weak shocks, 
is regarded as the most s uccessful energy source to b al- 
ance radiative cooling (see lMcNamara fc Nulsenll2007l for 
a recent review). Recent observations show that the 
AGN energies associated with X-ray cavities (~ 4PU, 
where P is the average pressure of the gas surrounding 
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the cavity, and V is the cavity's volume) are sufficient or 
nearly sufficient to stop the gas fro m cooling in most clus - 
ters containing detectable cavities dRaffertv et al.| [20061 . 
Furthermore, the mean power of these AGN outbursts 
increases in propor t ion to the cluster coo ling luminosity 
(jBirzan et al.l [200l iRaffertv et all |20M ) . This correla- 
tion suggests that AGN outbursts are triggered by the 
gas cooling and operates as a self-regulating feedback 
mechanism, which is shown to be essential in suppress- 
ing global thermal instability and thus in mai ntaining 
the ICM in the cool core state (|Guo et alj|2008f ). 

Despite these supportive clues, the detailed mecha- 
nisms that transfer the AGN jet energy to the ther- 
mal energy of the ICM are still far from clear. Nu- 
merical simulations have been extensively used to un- 
derstand this process. X-ray cavities are usually thought 
to form at the tips of AGN-induced jets. During cavity 
formation, weak shock wav es emerge as the gas is dis- 
placed, and heat the ICM (jBriiggen et all 120071) . The 
buoyantly-rising cavities may induce sound waves, which 
may be viscously dissipated in the ICM and thus heat 
the g as (|Ruszkowski et all 120041 : but see [Mathew s et ail 
2006). Radio synchrotron emission has been observed 
from many X-ray cavities (Birzan et al. 2004), suggest- 
ing that the cavities are formed by AGN-induced cosmic 
rays (CRs), which m ay leak into the ICM and heat the 
gas (IGuo fc Ohll2008l) . 

In most previous simulations, for simplicity, it is as- 
sumed that X-ray cavities are filled with ultrahot ther- 
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mal gas, but CR-filled cavities have received much less 
attention. The only explicit studies of the forma- 
tio n and evolution of X-ray cavities filled with CRs 
arelMathews fc Brighentl fl2008aI ). lMathews fe Brighentil 
(2008bl). and iMathewsl ()2009[ ). who show several inter- 
esting effects of CR-formed X-ray cavities: ICM heat- 
ing due to weak shock generation, ICM cooling asso- 
ciated with cluster expansion and outward mass trans- 
fer, and the dynamical relationship between post-cavity 
X-ray filaments and extended radio lobes. However, 
these studies focus on low-energy cavities (~ 10 58 — 10 59 
erg) in the Virgo cluster. In this paper, we apply the 
same model to study the remarkably large X-ray super- 
cavities observed in the cluster MS 0735.6+7421 (here- 
after MS0735) , created by the m ost en ergetic AGN out- 
burst known (IMcNamara et all l 2005t IGitti et al.l 120071 : 
IMcNamara et a l. 2009). In this cluster, two approxi- 
mately symmetric cavities are observed in X-ray images, 
each associated with relativistic gas emitting radio syn- 
chrotron emission. The energy released by this outburst 
is estimated to be 6.4 x 10 61 erg, assum ing an enthalpy 
of ApV per cavity (Raffcrt y et al J 1200 6). Weak shocks 
have also been detected in MS0735; a simple spherically- 
symmetric shock model suggests that the age and driving 
energy of the shoc k are t s = 104 Myr and E s = 5.7 x 10 6 
erg, respectively (|McNamara et all 120051). MS0 7 35 ha s 
been studied numerically by iSternberg fc Sokerl (2009), 
who assume that the cavities are formed by slow, mas- 
sive, wide jets with a total injected energy of 7.2 x 10 61 
erg and a very large injected mass ~ 1O 1O M0. 

One of the main objectives of the current paper is to 
investigate if AGN bubbles created by cosmic rays can 
reproduce the observed morphology of giant X-ray su- 
percavities in X-ray maps of MS0735. We show that the 
cavity morphology contains important information about 
the process of CR injection: CRs are likely to have been 
injected into the ICM at a series of locations (or contin- 
uously) along the direction of a je t as it moves outward, 
instead of at a fixed location (§ 13.11) . We then follow 
the time evolution of the ICM energies during and after 
cavity creation and study the energetics of X-ray super- 
cavities and their long-term influence on the cluster gas 
atmosphere. The investigation of this extremely powerful 
AGN outburst improves our understanding of how AGN 
outbursts affect the ICM, and their relevance in solving 
the cooling flow problem. 

As the largest virialized structures in the universe, 
clusters of galaxies are useful cosmological probes, pro- 
viding one of the current best constraints on the mean 
matter density O m , dark energy density J7de, and 
the dark energy equation of state parameter w (see 
lAllen et al.1 120081 and references therein). These con- 
straints are usually based on accurate measurements of 
the X-ray gas mass fraction profile f gas (r). In this paper, 
we study how the energetic AGN outburst in MS0735 af- 
fects the gas mass fraction by producing a buoyant out- 
flow and global expansion of the cluster thermal gas. X- 
ray supercavities also have a significant effect on the as- 
sumptions of hydrostatic equilibrium and spherical sym- 
metry, which are essential in measuring the cluster's total 
mass profile from X-ray observations. 

The rest of the paper is organized as follows. In § [51 
we describe basic time-dependent equations and our nu- 
merical setup. Our results are presented in § [3] We 



summarize our main results in § 2] with a discussion of 
the implications. 

2. EQUATIONS AND NUMERICAL SETUP 

To study the dynamical effects of CRs on the hot clus- 
ter gas following this extremely powerful AGN outburst 
in MS0735, we assume that AGN outbursts inject cosmic 
rays into the ICM, producing the obs erved radio lobes. 
X-ray supercavities and weak shocks (IMcNamara et al.1 
[20051 IGitti et al.l [20071 IMcNamara et al l 120091) . The 
combined evolution of the thermal ICM and cosmic rays 
may be described by the following four equations: 



dp 
~dt 



pV • v = 0, 



P^ = -V(P + P c ) - P V$, 



de 

— + V • (ev) = -PV • v, 



(1) 



(2) 



(3) 



de c 
dt 



V-(e c v) = -P c V-v + V-( K Ve c ) + S c , (4) 



where dj dt = d/dt + v • V is the Lagrangian time deriva- 
tive, P c is the CR pressure, e c is the CR energy density, 
K is the CR diffusion coefficient, S c is the CR source term 
due to the central AGN activity, and all other variables 
have their usual meanings. Pressures and energy densi- 
ties are related via P = (7 — l)e and P c = (7 C — l)e c , 
where we assume 7 = 5/3 and 7 C = 4/3. We do not in- 
clude radiative cooling since our intention is to study the 
evolution of X-ray supercavities in the cluster MS0735 
on timescales short compared to the multi-Gyr age of 
the cluster, i.e. we do not investigate the long-term bal- 
ance between radiative cooling and intermittent heating 
by AGN outbursts. The same set of equations and un- 
derlying p hysical assumptions ha v e been thoro ughly de- 
scribe d in iMathews fc Brighentil (|2008bf ) and IMathewsl 
( 2009) where further details can be found. Here we sim- 
ply summarize several modifications and reiterate some 
important points. 

Equations (fTJ) — (g| are solved in (r, z) cylindrical co- 
ordinates usin g a two-dimensional Eu lerian code similar 
to ZEUS 2D ([Stone fc Normanlll992f) : in particular, we 
have incorporated into the code a background gravita- 
tional potential, CR diffusion, and CR energy equation. 
The computational grid consists of 200 equally spaced 
zones in both coordinates out to 400 kpc plus an addi- 
tional 200 logarithmically-spaced zones out to 2 Mpc. In 
view of this large computational domain, we adopt re- 
flective boundary conditions for thermal gas and cosmic 
rays at the outer boundary as well as at the origin. 

For initial conditions, we adopt analytic fits to the de- 
projected electron nu mber density n e a nd temperature 
T profiles observed by IGitti et al J ([20071 ): 



n G (r) 



0.075 



0.01 



[l + (r/20) 2 ] 1 - 29 [l + (r/200) 2 ] 1 - 



cm 3 (5) 
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and 



T(r) = (Tf 1 ' 5 + T" 1 - 5 )" 2 / 3 keV, 



where 



Ti = 3.2 + 5.3(r/275) 



1.7 



T 2 = 8.5(r/275) 



-0.7 



(6) 
(7) 
(8) 

and r is measured in kpc. The fits are shown as solid 
lines in Figure 0] (a) and (b) . At the beginning of our 
simulation, the CR energy density is assumed to be zero 
throughout the cluster. The gravitational potential $ is 
set by assuming hydrostatic equilibrium at t = 0. 

X-ray cavities are usually thought to be inflated by 
bipolar jets emanating from an AGN in the central 
galaxy. The approximately symmetric double-cavity 
morphology in MS0735 suggests that the outburst in this 
cluster was nearly bipolar-symmetric, as we assume here. 
The jets deposit relativistic cosmic rays into small re- 
gions at their terminal points, which expand and form 
underdense bubbles producing the observed X-ray cav- 
ities. The injection of cosmic rays into the ICM is de- 
sc ribed in equation flU) by the s ource term S c . Simila r 
to IMathews fc Brighentil (|2008bD and IMathewsl (|2009D . 
we assume that the CRs are deposited into a Gaussian- 
shaped sphere of characteristic radius r s = 2 kpc located 
at r cav = (r,z) = (0,z cav ): 







; -[(r-r cav )/r„]^ 



when t < i agn , 
when t > i agn , 



0) 



where i agn is the duration of the CR injection (AGN ac- 
tive phase), and E agn is the total injected CR energy in 
one bubble (2E 3ign for the whole cluster). The integral 

of S c over space gives the the CR injection luminosity 
-EagnAagn associated with the creation of each bubble. 
To mimic the location and size of the observed X-ray 
cavity, we consider three cases for the location of the CR 
injection r cav : (a) z cav is fixed at 60 kpc (adopted in run 
MS-1); (b) z cav moves from 40 to 160 kpc at a constant 
speed during the CR injection phase (adopted in runs 
MS-2 and MS-2A); and (c) z cav moves from 40 to 160 
kpc with a constant deceleration during the CR injec- 
tion phase (the speed drops to zero at z cav = 160 kpc; 
adopted in run MS-3). See Table[T]for specific model pa- 
rameters in each run. The bulk kinetic energy of cosmic 
rays is negligible compared to the injected cosmic ray en- 
ergy E agn due to their negligible total rest mass. Here 
we assume that the energy that AGN jets deposit into 
the ICM is mainly in the form of cosmic rays, which are 
transported along the jets and created in strong shocks as 
the jets encounter the ICM. However, the energy content 
of real AGN jets is not clear and may also include mag- 
netic energy, the internal and kinetic energies of thermal 
(and superthermal) gas. An MHD jet simulation includ- 
ing cosmic ray acceleration (and re- acceleration) mecha- 
nisms is required to fully understand this process. 

In addition to their advection with the thermal gas, 
cosmic rays diffuse through the gas as described in 
equation (0}. The CR diffusion coefficient k is poorly 
known but may vary inversely with the gas density 
since the magnetic field is probably larger in denser gas 



TABLE 1 
List of Simulations. 







t b 

t'agn 


r 
1 cav 


''shock 


A f e 
ilJ shock 


Run 


(10 61 erg) 


(Myr) 




(10 s yr) 




MS-1 


4.5 


10 


a 


1.10 


1.41 


MS-2 


3.8 


10 


b 


1.16 


1.41 


MS-3 


3.5 


10 


c 


1.23 


1.41 


MS-2A 


3.8 


50 


b 


1.38 


1.34 



"The AGN energy released in the form of cosmic rays 
in one hemisphere during the outburst. 

6 The duration of the active AGN phase. 

c f C av = (0, 2cav) is the position where the cosmic rays 
are injected; see text in §[2] for three different ways of 
cosmic ray injection. 

d The time that the shock reaches r = 240 kpc along 
the semi-minor axis (i.e., the age of the energetic out- 
burst). 

e The Mach number of the shock along the semi-minor 
axis when it reaches r = 240 kpc. 



(IDolag et al.ll200l As in IMathews fe Brighentil (|2008bn 
and IMathewsl (|2009h . we adopt the following functional 
dependence of the diffusion coefficient on the gas den- 
sity: 



10 30 (n c0 /n o ) cm 2 s" 1 
10 30 cm 2 s" 1 



when n c > n c o, 
when n c < n e o, 



(10) 



where we take n o = 10 5 cm 3 in the rest of the paper. 
M odels with different va l ues of n r n have been explored 
bv lMathews fe Brighentil (|2008bD and the results are usu- 
ally not very sensitive to it. During their diffusion, cos- 
mic rays interact with magnetic irregularities and Alfven 
waves, exerting CR pressure gradients on the thermal 
gas (equation [2]). The shock generation and early evo- 
lution of cavities are mainly caused by the expansion of 
the surrounding gas driven by the CR pressure unless the 
coefficient of CR diffusion is so large (e.g., n o > 10~ 2 
cm~ 3 ) that it significantly affects the morphology of cav- 
ities. We neglect other more complicated interactions of 
cosmic rays with thermal gas, e.g., Coulomb interactions, 
hadronic collisions, and hydromagnetic-wave-mediated 
CR heating, which depend on the cosmic ray energy spec- 
trum and provide additional heating effects for the ICM 
(e.g., IGuo fe Oh|[200l . We defer a thorough study of 
these effects to future work. 

3. RESULTS 

3.1. Morphology and Structure of X-ray Supercavities 

Chandra observations show that the cluster MS0735 
hosts two giant cavities. Each cavity is slightly elon- 
gated in the radial directi on ('prolate') and is roughly 20 
kpc in (radial) diameter (jMcNamara et al.l [20051 120091) . 
Weak shocks have also been detected; the shock front 
is roughly elliptical with a semimajor axis of ~ 360 kp c 
and a semiminor axis of - 240 kpc (|Gitti et al.l [20071. 
Along the semiminor axis (perpendicular to the cavity 
axis), the shock Mach n umber is estimated to be ~ 1.41 
(jMcNamara et al.ll2005l ). 

One of the primary goals of this study is to investigate 
whether AGN bubbles created by cosmic rays can repro- 
duce these observational features. Adopting the radiative 
coolin g function n c n{ A(T, Z) from lSutherland fc Dopital 
(1993), we project the emissivity and produce synthetic 
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Fig. 1. — Synthetic X-ray surface brightness maps (projections of n c niA(T, Z) viewed perpendicular to the bubble symmetry axis) for runs 
MS-1 (top), MS-2 (middle), and MS-3 (bottom) in cgs units with logarithmic scale. Left column: t = 50 Myr; right colu mn: t = t s h oc k (see 
Table [ll. X-ray supercavities are clearly visible. The prolate shapes of the cavities in runs MS-2 and MS-3 fit observations (McNamara ct al. 
2009) better than the oblate cavity seen in run MS-1. 

dially prolate cavities, consistent with observations. Our 
methods of CR injection are probably idealized, but the 
resulting cavity shapes strongly suggest that CRs may 
indeed be injected into the ICM at a scries of locations 
(or continuously) as AGN jets move outward, instead of 
at a fixed location as assumed by most au thors (e.g., 
IRuszkowski et al.ll2007t IBriiggen et al.M2009f >. 

Images of the cluster MS0735 combining both X-ray 
and radio wavelengths indicate that the X-ray cavities 
are filled with r a dio e mission, as seen in Fig. 1 of 
iMcNamara et all (|2009f) . While the overall radio im- 
age in MS0735 has an approximate north-south bipolar 
symmetry relative to the cluster-centered galaxy, several 
deviations from perfect symmetry are apparent. Weak 
radio emission observed within about 50-70 kpc from the 
center has a noticeably different alignment. It is possi- 
ble that the inner radio structure has been moved from 
an original north-south alignment by flows in the hot 
gas or that MS0735 has ejected radio plasma in several 
different directions. Most of the bright extended radio 
emission comes from radially elongated elliptical-shaped 
regions that are well separated from the cluster center. 
Considering a typical synchrotron-loss timescale of less 
than about 100 Myr (depending on both the magnetic 
field strength and electron Lorentz factor) for relativistic 
electrons in clusters, the decline in radio flux close to the 
center of MS0735 could arise if the radio electrons de- 



X-ray maps of our runs at different times. Here n; is 
the total number density of ions and we ta ke an average 
metallicity of Z = 0.4 (see lGitti et al.ll2007t) . m is related 
to the proton number density n p via fij = l.ln p , and 
thus the molecular weight is /i = 0.61. Figure ([1]) shows 
synthetic X-ray surface brightness maps in logarithmic 
scale for our three typical models (MS-1, MS-2, MS-3) 
as listed in Tabled] at two different times t — 50 Myr (left 
panels) and t — t s i loc k (right panels), where i s hock is the 
time when the shock propagates to r = 240 kpc along 
the semiminor axis, determining the age and energy of 
the outburst. X-ray deficient cavities and shock fronts 
are clearly seen in all these runs, which differ mainly in 
the location of cosmic ray injection. The CR injection 
time in all these runs is assumed to be 10 Myr, and the 
injected CR energy E agn is determined so that the shock 
along the semiminor axis has a Mach number of 1.41 
when propagating to r = 240 kpc. 

In run MS— 1, cosmic rays are injected at a fixed posi- 
tion (z cav = 60 kpc in the z axis). The top panels of Fig- 
ure (flj show that the X-ray cavity has a roughly oblate 
shape. However, the observed cavity is elongated along 
the radial direction, which suggests that CRs may be in- 
jected into the ICM at a range of z cav as the jet tip moves 
outward. Thus, in runs MS-2 and MS-3, we assume that 
z cav moves from 40 to 160 kpc within t < t agn (see §[2]for 
details) . The middle and bottom panels clearly show ra- 
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posited there have radiated most of their initial energy. 
This is consistent with our result that CRs are injected 
into the ICM at a series of locations (or continuously) 
as AGN jets move outward. Since electrons with lower 
energies have longer synchrotron-loss times and emit in 
longer wavelengths, future radio observations at different 
wavelengths could provide further insights regarding the 
formation and evolution of the MS0735 X-ray cavities. 

Slices of logarithmic gas density through the cluster 
center are shown in Figure [2] for runs MS-1 (top pan- 
els) and MS-2 (bottom panels) at t — 50 Myr (left) and 
t = tshock (right). Contours show the logarithmic CR en- 
ergy density. The bubbles are clearly seen as low-density 
regions along z axis, where cosmic rays also concentrate. 
In run MS-2, where CRs are injected from z cav = 40 to 
160 kpc along z axis, the bubble structure at t = i s hock is 
much more complex than that in run MS-1, although the 
cavities in the projected X-ray surface brightness maps 
are simple and smooth in both cases. Since the prolate 
cavity in run MS-2 fits observation better, this raises the 
possibility that the cavities in X-ray maps may have a 
complex structure, with large variations of gas density 
and CR energy density at different locations inside cavi- 
ties. Also of interest is the buoyant flow of the CR cav- 
ities. At the time of observation i s hock the dark cavities 
in Figure 2 have moved well away from the CR source 
regions, z/2 = 30 kpc for MS-1 and 20 < z/2 < 80 kpc 
for MS-2. 

In addition to the shock confirmed by observation, an- 
other shock can be seen in each panel of Figure Q] and on 
closer inspection also in Figure [5] This second shock, ob- 
viously due to the reflection of the first shock at the sym- 
metry boundary z = 0, represents the shock produced by 
the other identical cavity in the cluster at z < 0. It goes 
through the already shocked gas and has a much smaller 
Mach number (with a maximum value of ~ 1.1 far from 
the X-ray cavity). Observational efforts to detect this 
second shock (even non-detection) could provide inter- 
esting information on the geometry and simultaneity of 
bipolar jets and cavities in clusters. 

3.2. Outburst Momentum and Black Hole Recoil 

As the most powerful known AGN outburst, MS0735 
highlights another fundamental concern regarding the 
joint evolution of massive elliptical galaxies and their 
central supermassive black holes - why the masses of 
the black holes stellar bul ges are so strongly co rrelated 
(jFerrarese fc Merrittil2000t ' iGebhardt et al.ll200Ct ). Since 



energy M Ci 



3.8 x 10 ergs contained in each large 



cavity of MS0735 is purely relativistic with a mass equiv- 
alent of M cav = £ ag „/c 2 = 2.1 x 10 7 M Q . If this energy 
was supplied to the cavity by a relativistic jet as obser- 



elliptical galaxies are commonly regarded as the end 
products of galaxy mergers, the black hole-stellar mass 
correlation suggests that the black holes merge in the 
same fashion as the surrounding stars. However, de- 
tailed studies show that asymmetric gravitational radia- 
tion emitted during the final stages of black hole merg- 
ers results in recoil blac k hole velocities < 4000 km s _1 
(|Campanelli et all 120071) sufficient to displace th em far 
from the galactic cores (Madau & Quatacrt 2004). Nev- 
ertheless, dedicated searches for peculiar AGN recoil ve- 
locities or binary AGNs have resulted in surprisingly few 
candidates. 

However, active, jet-emitting black holes in powerful 
radio galaxies like MS0735 could also receive very large 
recoil kicks if the jets are not exactly balanced by an 
identical counterjet. Suppose for example that the total 



vations indicate, its momentum p = M c , 



1.26 x 10 51 



gm cm s 1 is sufficient to produce a recoil velocity of 
6300 km s" 1 for a black hole of mass M bh = 10 9 M . 
This velocity is large enough to completely eject the mas- 
sive black hole fr om typical group an d cluster-centered 
elliptical galaxies (|Merritt et al.l l2004). To avoid this un- 
desired fate, it is necessary that the momentum of the 
MS0735 jets, and those of other powerful radio galaxies, 
are created in almost identical oppositely directed jet and 
counter-jet pairs during the entire period of intense jet 
activity which may last tens of million years. 

While a pair of identical opposing jets may not al- 
ways result in identical cavities, it is not unreasonable 
to expect such cavity pairs to be quite common, at least 
in their early development. If so, we would expect the 
shocks from each young cavity to collide at the symmetry 
plane perpendicular to the jet axis, reflecting as an addi- 
tional outgoing shock similar to those shown in Figure ([T]) 
which should be visible in deep X-ray images. However, 
X-ray images of nearby clusters like Perseus and Virgo - 
as well as more distant clusters like MS0735 - often show 
cavities or cavity-related thermal features that appear to 
have no counterpart or misaligned counterparts in the 
opposite hemisphere. 

In particular, the image of MS0735 in lMcNamara et al.1 
(2009) shows that the two radio lobes and their associ- 
ated X-ray cavities are not mirror images of each other. 
The most plausible explanation for this is that the cavi- 
ties have been pushed aside and distorted by motions in 
the nearby hot cluster gas. We estimate from this figure 
that the cavities have been displaced transverse to their 
common jet axis by at least a distance of d ~ 70 kpc. 
At the location of the radio cavities the sound speed 

in the hot gas is c s = (jk^T/ [im v ) x ^ 2 = 513T^ e y km 
s- 1 « 1150 km s- 1 (for T koV = k h T/l keV « 5) and the 
hot gas density is 9.5 x 10~ 27 gm cm" 3 . Can the non- 
axisymmetric radio emission dislocations seen in MS0735 
be generated by subsonic flows in the ambient gas? To 
create the observed spatial perturbations in MS0735, 
consider a ram pressure pi> 2 am acting on the radio lobe of 
area A = nab (where a ~ 115 kpc and b ~ 70 kpc for the 
northern lobe in MS0735) which accelerates the lobe by 
2d/i 2 am w pu 2 am A/M cav in time i ram sa i shock « 10 8 yrs. 
The required gas velocity, u ram ~ 9 km s _1 is very much 
less than the sound speed, suggesting that purely rela- 
tivistic radio lobes are extremely vulnerable to distortion 
by very modest flows in the ambient cluster gas. 

However, it is likely that the jet that produced the 
radio emitting cavities in MS0735 acquired significant 
masses of very hot but non-relativistic gas during their 
passage to the cavity. Suppose for simplicity that the 
pressure of this non-relativstic gas dominates in the 
cavities and that the cavities are in pressure equilib- 
rium with the local gas, p C av?cav = pT. The ellip- 
soidal cavity mass is now M cav = p cav (4/3)7ra6 2 and the 
ram velocity required to move the cavity a distance d, 
<W « [(8/3)(T/T cav )6dA 2 am ] 1 /2 ; i s 350 and 110 km s' 1 
for T cav = 50 and 500 keV respectively. As before, it is 
plausible that the cluster gas in MS0735 has subsonic ve- 
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Fig. 2. — Slices of log (p/g cm -3 ) in runs MS-1 (top) and MS-2 (bottom). Left column: t = 50 Myr; right column: t 

— ^shock (see Table 

[TJ. White contours showing log (e cr /erg cm -3 ) indicate that cosmic rays concentrate within the X-ray supercavities. 



locities of this magnitude, so the lack of mirror symmetry 
in the radio sources tells little about their mass content 
or the almost perfect cancellation of the jet momentum 
that created them. Nevertheless, a statistical study of 
the distortion and transverse displacement of radio lobes 
or X-ray cavities would contain information about clus- 
ter gas velocity fields and the baryonic component in the 
radio cavities. 

3.3. Energetics and Heating Effects 

Since the morphologies of the cavity and shock front 
in run MS-2 agree with observations very well, we take it 
as our fiducial model to study the effects of the energetic 
AGN outburst on the cluster MS0735 in the rest of the 
paper. 

The top panel of Figure [3] shows time evolution of 
global energies integrated over the whole cluster in run 
MS-2. Ei n j(t) (solid lines) represents the total injected 
CR energy until time t, which equals to the overall in- 
crease of cluster energies, including the gas kinetic energy 
Ek (short dashed lines), the change in gas thermal energy 
Ai?th (dotted lines), the change in gas potential energy 
A-Epot (dot dashed), and the total CR energy E cr (long 
dashed lines). Clearly, E CT (t) is much less than the to- 
tal injected CR energy E-^(t) during the CR injection 
phase, indicating that CRs lose a significant fraction of 
their initial energy due to pdV work and shock genera- 
tion. At time t = i agn (the left dotted vertical line) , more 



than 60% of the injected CR energy has been converted 
into thermal and kinetic energies as the thermal gas is 
shocked. The CR energy lost during the injection phase 
is expected to be less if the same amount of energy is 
injected over a longer timescale. To check this explicitly, 
we performed another run MS-2A for comparison, which 
is the same as run MS-2, except that the CR injection 
timescale is much longer (t a gn = 50 Myr). The resulting 
energy evolution, plotted in the lower panel of Figure [31 
confirms that less (~ 50%) of the injected CR energy has 
been converted into gas thermal and kinetic energies at 
the end of the injection phase. We therefore conclude 
that a large fraction (> 50%) of the AGN energy is de- 
livered to the ICM during the formation of X-ray cavities 
through shock generation driven by pdV work. This frac- 
tion increases with the CR luminosity E agn /t agn during 
the CR injection phase. 

The evolution of the gas thermal and potential ener- 
gies in Figure |3j suggests that the thermal energy of the 
ICM responds to cavities in two competing ways. On the 
one hand, cavities directly heat the ICM through weak 
shocks, as clearly shown by the increase of thermal en- 
ergy with time during t < 100 Myr in Figure [3J On 
the other hand, as cosmic rays displace cluster gas, the 
entire cluster expands (see Fig. [5] and § 13.41 for more 
details) and thus cools as the gravitational potential en- 
ergy of the cluster increases. Local heating dominates 
during the early stage of cavity evolution, while expan- 
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Fig. 3.— Global energy evolution in runs MS-2 and MS-2A. The 
energies are calculated in one hemisphere and are labeled as fol- 
lows: the injected cosmic ray energy E in j (solid), change in thermal 
energy AE til (dotted), cosmic ray energy E CI (long dashed), kinetic 
energy (short dashed), and change in potential energy AE po t 
(dot dashed). The vertical lines in each panel show the time t ag n 
when the CR injection ends (at left) and ts^ocfe when the shock 
propagates to r = 240 kpc (at right), respectively. 

sion cooling dominates during the later stage when the 
shock becomes much weaker. This is indicated by the 
decline of A-Bth and the accompanying increase of grav- 
itational energy during t > 100 Myr in Figure |3] This 
'global cooling' effect has be en explicitly discussed by 
lMathews""fe Brighcnti (2008b) for a much less powerful 
AGN outburst in Virgo cluster (10 58 erg); our new re- 
sults for MS0735 suggest that it is likely to be a universal 
feature of AGN outbursts regardless of their energy. 

Figure [4] shows radial profiles of emission-weighted 
spherically-averaged gas quantities, including (a) elec- 
tron number density, (b) temperature, (c) pressure, and 

(d) gas entropy (S = k B T/nl /3 ) in runs MS-2 and MS- 
2A. At t — ishock, the gas density in inner regions, which 
are already affected by the AGN activity, drops in both 
models, confirming that the cluster is undergoing global 
expansion. Consequently, the expansion of these regions 
results in a drop in the gas temperature, as seen in the 
short dashed line in Figure [4Jd for run MS-2A. In run 
MS-2, the CRs are injected in a much shorter timescale, 
resulting in stronger shock heating, which surpasses the 
effect of global expansion at time i s hock- Thus the central 
cluster gas temperature and entropy both increase in run 
MS-2 (see Fig. gU). At later times, t > 2t shock « 20t agn , 
cooling by global expansion begins to do minate. As ex- 
plained by Mathews fc Brighentil (|2008bf ) , the expansion 
of cluster gas lasts much longer than the cavity since the 
volume V of gas displaced by the initial cavity is approx- 
imately preserved after the CRs diffuse widely into the 
cluster gas. 

3.4. Effects on the Gas Mass Fraction 

By assuming that clusters contain a proportional share 
of all cosmic baryons, cluster gas mass fractions (/ gas ) 
have been used as a probe of the universal ratio of 
baryon to total matter densities, fib/^m, providing one 
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Fig. 4. — Radial profiles of emission-weighted spherically- 
averaged (a) electron number density ra e , (b) temperature T, (c) 

pressure P, and (d) entropy S (S = k-gT/n^ 3 ) in runs MS-2 and 
MS-2A. The long-dashed and dot-dashed lines in panel (c) repre- 
sent the cosmic ray pressure distribution in run MS-2 and MS-2A, 
respectively. 



of the current best constraints of f2 m when combined with 
accurate determinations of fib from cosmic nucleosyn- 
thesis calculations or observations of cosmic microwave 
background (CMB) anisotropics. Furthermore, measure- 
ments of the apparent evolution of / gas with rcdshift 
can be used t o constrain the da rk energy content of the 
universe (see lAllen et al.l I2008I and references therein). 
While cluster baryon mass fractions seem to provide 
convincing cosmological constraints, systematic uncer- 
tainties, especially the apparent effects of cooling and 
heating on the distribution of baryons within clusters, 
must be well understood and controlled (see relevant 
discussions in iBorgani fc Kravtsovl [2009) . It has been 
noticed that cluster baryon fractions estimated from X- 
ray observations are significantly lower than the universal 
baryon fraction inferred from WMAP data (|Ettorill2003t 
iVikhlinin eHa! I2006I : [McCarthy et al.l 120071 ): in partic- 
ular. IVikhlinin et al.l (|2006l ) found that for a sample of 
13 Chandra clusters, gas mass fractions within r25oo are 
around 0.05-0.10, a fa ctor of 1.7-3.3 below the universal 
baryon fraction 0.167 (|Komatsu et al.ll2009ft . Such a dis- 
crepancy cannot b e contained in the stel lar /cool baryon 
contribution (e.g.. iMcCarthv et al.ll2007j ). 

However, the baryon fraction profile /bar(^) can be cre- 
ated and shaped by AGN outbursts. Figure [5] shows 
the time evolution of the cumulative spherical gas mass 
M gas (r) (top panel) and of its fractional change (bot- 
tom panel) in run MS-2, where M gas (r) is defined as the 
total gas mass within a cluster-centric radius r. Compar- 
ing gas mass profiles at different times, one can clearly 
see that, as a result of cavity formation, the cluster un- 
dergoes a global expansion and gas mass is gradually 
transported outward. After t ~ 300 Myr, the radial 
profile of M gas (r) in the central regions (within ~ 200 
kpc) becomes quite steady and within ~ 100 kpc about 
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20% — 30% of the initial gas mass has been transported to 
outer cluster regions. At t = 1 Gyr, around 4% of the gas 
mass originally within 500 kpc, ~ 6 x 10 11 M©, has been 
transport ed outward. No te that, r25oo ~ 465 kpc for 
MS0735 (|Gitti et alj|2007h . During the lifetime of clus- 
ters, a series of intermittent AGN outbursts are expected 
to be triggered, preventing the formation of strong cool- 
ing flows and transporting a significant fraction of the 
cluster gas to large radii. This large post-cavity outflow 
can explain most of the discrepancy between the cluster 
baryon mass fraction and the universal baryon fraction 
inferred from WMAP data. Without correcting for this 
effect, systematic errors inevitably occur when using / gas 
at r 2 5oo ( ^25no is typically ~ 300-60 kpc for galaxy c lus- 
ters; see IVikhlinin et all [20061 and I Allen et all 120081 ) to 
constrain cosmological parameters; for instance, Q m will 
be systematically overestimated. 

An interesting feature of Figure [5] is that the post- 
cavity fractional change of the gas mass decreases with 
radius, well beyond the 40 < r < 160 kpc CR source re- 
gion for run MS-2. This directly results in the trend that 
the gas mass fraction / gas (< f) within a cluster-centric 
radius r increases with r, which is consis tent with obser- 
vations (e.g., Fig. 7 of lPratt et al.1120091 ). This may also 
naturally explain the observational fact that the cluster 
gas mass fraction within ^2500 (or 7-500) increases with 
clust er temperature (e.g., see Fig. 21 of IVikhlinin et all 
2006). For higher-temperature clusters, r 2 5oo (or 7-500) is 
usually larger, and thus a smaller fraction of the baryon 
gas is transported beyond this radius if the strength of 
AGN outbursts does not increase propo rtionally to over- 
take t his effect. On the other hand, iMcCarthv et al.l 
(2007) suggested that the discrepancy between the clus- 
ter baryon mass fraction and the universal baryon frac- 
tion inferred from WMAP data is mainly due to the fact 
that the actual value of f2 m is higher than the best-fitting 
WMAP value. We argue that this explanation would re- 
sult in a constant gas mass fraction across all clusters, 
which is apparently inconsistent with observations. For 
galaxy groups, observations show that the gas mass frac- 
tion at T25oo is around 0.05, considerably smaller than 
the typical value of ~ 0.09 for clusters. The gas frac- 
tion als o increases with the mean temperature in galaxy 
groups (ISun et all 12009ft. 

It was arg ued bv IMcCarthv et al.l (|2007f ) that the total 
energy required to produce this discrepancy by trans- 
porting gas mass outward may be extremely large (~ 
10 63 erg). However, their simple calculations assume ar- 
bitrary extrapolations of gas density profiles to cluster 
outer regions not yet observed, where the transported 
gas, which determines the required increase in the gas 
potential energy, is located. On the other hand, our sim- 
ulations directly show that single strong AGN outbursts 
transport a significant amount of gas outward, and the 
increase of potential energy during the cluster expansion 
is accompanied by the decrease of thermal energy, as gas 
moves outward down the thermal pressure gradient (see 
the evolution of thermal and potential energies in Fig. 3). 
Consequently, the total AGN energy needed to explain 
the gas mass fraction discrepancy, i.e., the net increase 
of the total cluster energy, may not be energetically pro- 
hibitive. Further investigation of this is obviously beyond 
the scope of the current paper, and is thus left to future 
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Fig. 5. — Time evolution of the cumulative spherical gas mass 
within radius r (top panel) and of its fractional change (bottom 
panel) in run MS-2. Here M gaSi o(r) stands for Mg as (r) at t = 
Myr. After t ~ 300 Myr, the radial profile of Af gas (r) in the central 
regions (within ~ 100 — 200 kpc) becomes quite steady and around 
20% — 30% of the gas mass has been transported to large radii. 



work. 

The effect of AGN outbursts on the gas fraction has 
been studi ed with cosmological simulations of galaxy 
clusters by lPuchwein et al. (2008). Using a simple model 
for AGN mechanical heating, they find that AGN heating 
significantly reduces gas fraction within r§oo in groups 
and poor clusters while having much less er effects on 
massive clusters. The analytical study by iPopej (2009) 
also suggests that AGN energy injection redistributes 
the ICM within the cluster central regions. These find- 
ings are consistent with our detailed simulations: the ef- 
fect of AGN outbursts on gas fraction decreases with ra- 
dius, and becomes quite small a t rspp for massive clusters 
(r 500 ~ 1340 kpc for MS0735: IGitti et all 12007ft . How- 
ever, observations indicate that / ga , s (< r) incr eases with 
r at least up to r ~ r 50 o (jPratt et al.ll2009t ). which is 
roughly the largest radius for which Chandra and XMM- 
Newton data require no model extrapolation. It may sug- 
gest that the central AGN produces some X-ray cavities 
at radii much larger than what we assume here, which 
transport a large amount of gas out of rsoo, and are dif- 
ficult to detect due to low X-ray surface brightness in 
these outer regions. Alternatively, it is also possible that 
massive clusters are formed from mergers of less massive 
systems, which have already lost a significant amount of 
gas by AGN outbursts. Further studies are required to 
determine which scenario dominates. 

3.5. Assumptions of Hydrostatic Equilibrium and 
Spherical Symmetry 

By assuming hydrostatic equilibrium and spherical 
symmetry, radial profiles of cluster gas density and tem- 
perature have been widely used to measure the to- 
tal gravitational mass within radius r, M tot (r) (e.g., 
IGitti et alJl2007ft : 

k h Tr f d\nn c d\nT \ 
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Fig. 6. — Radial profiles of the total gravitational mass Mtot(r) 
in run MS-2 calculated at several times by assuming hydrostatic 
equilibrium and spherical symmetry. The dashed and solid lines 
largely overlap. 

where G is the gravitational constant, fcb is the Boltz- 
mann constant, m p is the proton mass, and /i = 0.61 
is the molecular weight. Although this method provides 
one of the current best measurements of mass profiles for 
relaxed clusters, the underlying assumptions, i.e., spheri- 
cal symmetry and hydrostatic equilibrium, must be care- 
fully tested. Systematic errors, in particular the non- 
thermal pressure support from CRs or turbulent motions, 
have been ex tensively investiga ted by cosmological sim- 
ulations (see Nag ai et al.l 120071 and references therein). 
Here we study how X-ray supercavities, which have yet to 
be included in cosmological simulations, affect the mea- 
surements of cluster mass profiles found using equation 

OH)- 

For simplicity, we ignore self-gravity of the thermal 
gas. Thus the total gravitational mass, which is domi- 
nated by dark matter, does not evolve with time in our 
simulations. Figure |5] shows radial profiles of the to- 
tal gravitational mass M tot {r) in run MS-2 calculated 
at various times using equation (|11|) and the emission- 
weighted spherically-averaged gas density and temper- 
ature profiles (see Fig. SJ. The solid line in Figure |6] 
represents the cluster mass profile before the AGN out- 
burst, while the dotted line shows the mass profile at 
t = ishock = 1.16 x 10 8 yrs when the shock propagates 
to r = 240 kpc along the semiminor axis. By comparing 
these two lines, one can clearly see that the total clus- 
ter mass measured by equation (|11|) at t = i s hock is not 
very accurate for r ~ 30 — 300 kpc, which is exactly the 
region most strongly affected by the X-ray supercavities 
(see the right-bottom panel of Fig. [2]). As the cavity 
rises in the cluster atmosphere, the ICM near the cavity 
is significantly disturbed, which can be easily seen in the 
spherically-averaged radial temperature and density pro- 
files (the top two panels of Figure 2]). Thus, both the as- 
sumptions of hydrostatic equilibrium and spherical sym- 
metry are violated, and equation (jlip does not hold very 
well in the cavity region when t ~ i s hock- Since t — £ s hock 
is also the time when the deprojecte d gas temperatur e 
and density profiles were observed bv lGitti et al.l (|2007l ). 
the varying slope dM to t/dr visible in our initial mass 
distribution (solid line in Figure [5]) near the cavity re- 
gion also deviates somewhat from the expected N avarro- 
Frenk- White (NFW) form (jNavarro et all 119971) . This 



unavoidable deviation does not affect our conclusions 
which are all based on the differential changes in the clus- 
ter gas due to the cavities. The dashed line in Figure [6] 
represents the cluster mass measured by equation (JTTJ) at 
a much later time, t = 1 Gyr, when the X-ray supercav- 
ity has already been destroyed. At this time, the dashed 
line agrees very well with the solid line in regions r < 700 
kpc, indicating that the cluster has returned to hydro- 
static equilibrium and spherical symmetry and equation 
(jllj) provides a very accurate measure of the cluster mass. 
At r ~ 1000 kpc in Figure [6] the dashed line oscillates 
around the solid line due to outward-propagating weak 
shocks produced by the outburst. Evidently, more accu- 
rate total mass determinations for MS0735 using equa- 
tion (1111) could be found from observations of n e (r) and 
T(r) in cluster quadrants that avoid the cavity region. 

In summary, our simulations show that X-ray cavities 
produced by AGN outbursts may significantly disturb 
the cluster gas near the cavities and the assumptions 
of hydrostatic equilibrium and spherical symmetry may 
thus not hold very well in these regions. This results 
in systematic errors in measuring the cluster mass using 
spherically-averaged profiles of gas density and temper- 
ature. For the powerful outburst in the cluster MS0735, 
the typical error in -Mt o t(0 at radii near the cavities is 
around 10 — 30%, but may be much larger in some lo- 
cal regions. We have also done similar calculations with 
much less energetic AGN outbursts, and found that the 
cavities created are smaller and thus affect smaller cluster 
regions. X-ray cavities do not affect the measurements 
of cluster mass in the regions far away from them, in 
particular, in the cluster outer regions (e.g., rsoo). 

4. SUMMARY AND DISCUSSION 

By conducting a suite of two-dimensional axisymmct- 
ric hydrodynamical simulations, we investigated the for- 
mation and evolution of the X-ray supercavities recently 
observed in the cluster MS0735, which is the most ener- 
getic AGN outburst known. We assume that the X-ray 
supercavities are inflated by cosmic rays injected from 
AGN jets, and follow the co-evolution of the cluster gas 
and the CRs. We then study the thermal and hydro- 
dynamic effects of this energetic AGN outburst on the 
cluster MS0735. Here we briefly summarize our main 
findings: 

1. X-ray deficient cavities and weak shocks are success- 
fully created by AGN-generated cosmic rays. In partic- 
ular, run MS-2 reproduces very well many characteristic 
observational features of the cluster MS0735, such as the 
size and location of the cavities, the location and strength 
of the shock. In this run, the age and total energy of this 
AGN outburst are 116 M yr and 7.6 x 10 61 erg, consiste nt 
with previous estimates (jMcNamara et al.M2005l I2009T >. 

2. Assuming that X-ray cavities are generated by rel- 
ativistic AGN jets, the jet that produced the north cav- 
ity in run MS-2 had a momentum ~ 2 x 10 7 M Q c which 
must be balanced by an almost identical southerly jet 
ejection. Otherwise, a one-sided jet in MS0735 would 
cause its AGN black hole to attain a recoil velocity of 
~ 6000(Mbh/10 9 M Q ) _1 km s" 1 , which is large enough 
to completely eject a massive black hole from the cen- 
ter of even the most massive elliptical galaxies and 
their surrounding g roup or cluster dark matter halos 
(|Merritt et al.l 120041) . The radio emission from the cav- 
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ities in MS0735, while clearly bipolar, is not perfectly 
symmetric. This could be due to the production of some- 
what different cavities, perhaps at different cluster radii, 
by two jets that were identical when they left the AGN. 
In addition, we show that small subsonic motions in the 
hot gas of MS0735 are sufficient to explain the axisym- 
metric cavity radio and X-ray emission observed. How- 
ever, if the cavities formed by the two opposing jets were 
initially approximately equal in power and formed at sim- 
ilar distances from the cluster center, the shocks expand- 
ing away from each cavity would collide at the symmetry 
plane, producing a second, somewhat weaker shock in 
each hemisphere as seen in our Figure [T] Efforts should 
be made to detect these symmetry shocks in deep X-ray 
images since they contain information about the identical 
nature of the two jets and their cavities. 

3. Run MS-1, where CRs are injected at a fixed loca- 
tion, produces an oblate, radially-flattened X-ray cavity 
unlike the observed cavity which is radially elongated. In 
contrast, in run MS-2 which fits observations very well, 
the CRs are injected continuously along the jet direction 
as it moves away from the cluster center. The radial elon- 
gation of the cavity may also be produced if the CRs arc 
injected at a series of locations along the jet direction. 

4. During the CR injection phase, more than half of the 
injected CR energy is converted into thermal and kinetic 
energy of the ICM. The CRs lose energy thorough pdV 
work and shock generation as they displace the cluster 
gas. For instance, ~ 63% of the injected CR energy has 
been lost at t = i agn = 10 7 yrs in run MS-2. The fraction 
of initial CR energy lost at time £ agn increases with the 
CR luminosity during the injection phase. 

5. The thermal energy of the ICM responds to CR- 
generated cavities in several complex ways. As the cav- 
ities are inflated in MS0735, shocks are driven into the 
surrounding gas, heating the local ICM. Mixing by CR 
buoyancy may also increase the gas entropy near the 
post-cavity site as relatively low entropy, CR-enriched 
gas near the cluster core flows outward with the CRs 
and is replaced by inflowing gas of higher-entropy. On 
the other hand, the production of new cavities causes the 
cluster gas to readjust outwards in the cluster, increas- 
ing its global potential energy and decreasing its thermal 
energy (cooling) as its density decreases. Local heat- 
ing by cavity shocks dominates during the early stage 
of cavity evolution, whereas cooling by global expansion 
dominates at later times, when the total cluster ther- 
mal energy decreases below the original cluster thermal 
energy before the outburst. Thus, from a global perspec- 
tive, the cluster gas is eventually cooled by the cav ity, as 
previously shown bv lMathews fc Brighentil (|2008b[ ) for a 
low-energy outburst in the Virgo cluster. We have now 
confirmed this for the cluster MS0735 which is undergo- 
ing the most powerful AGN outburst known. This ulti- 
mate allocation of cavity energy is likely to be universal 
for all AGN outbursts. 

6. The creation of X-ray supercavities in MS0735 with 
CRs produces a huge, long-lasting outward migration 
of cluster gas to large radii. This mass outflow arises 
because of the global cluster expansion due to the gas 
displaced as the cavity forms and heating by outward- 
propagating shock waves driven by the expanding cavity. 
The mass outflow following the single energetic outburst 
in MS073 is large. Repeated AGN outbursts may nat- 



urally explain why the cluster baryon mass fraction in 
MS0735 and other similar clusters is much less than the 
universal baryon fraction inferred from WMAP data, and 
why it increases with the cluster temperature. 

7. X-ray supercavities significantly disturb the cluster 
gas near them, locally invalidating the assumptions of 
hydrostatic equilibrium and spherical symmetry. This 
results in systematic errors (typically around 10 — 30% for 
MS0735) in measuring the cluster mass in these regions 
using spherically-averaged gas density and temperature 
profiles. 

8. Assuming that X-ray cavities are filled with rel- 
ativistic CRs, APV is often used to estimate the en- 
ergy released by each AGN jet. When combined with 
a typical (buoyant) cavity lifetime t\a e , the jet power 
APV/ tiife can be estimated and compared with the bolo- 
metric X-ray luminosity of the cluster gas (jBirzan et al.l 
l200l . Since the total CR energy contained inside the 
cavity is E CT V = 3PV, where E a is the average CR en- 
ergy density in the cavity, about 1/4 of the jet energy 
(APV - E a V = PV) is delivered to the ICM during the 
cavity formation in these approximations. In our simula- 
tions, we directly studied the energetics of the cavity, and 
find that more than half of the injected CR energy goes 
into the ICM, i.e., the CR energy lost during the cavity 
formation is actually more than 3PV. This fraction in- 
creases with the CR luminosity. Although APV provides 
an approximate estimation for AGN jet energy, the ac- 
tual AGN energy may be a few times larger (typically 
~ 6PV — 10 PV). This is consistent with previous simu- 
latio ns of cavities produc ed by (non-relativistic) thermal 
jets ()Binnev et al.l l200~7h . which, due to the additional 
inclusion of a large amount of kinetic energy in thermal 
jets, find that AGN energy is around 15 PV, even larger 
than our estimation. Our calculations further show that 
APV begins to decrease after the c avity formation as CRs 
diffuse through the cavity walls (|Mathews fc Brighentil 
2008b). More generally, if APV/t\\f e underestimates the 
cluster jet power, it becomes easier to match this heating 
with the bolometri c rad iative losses Lx as p roposed by 
IBirzan et all (pOOl and lRaffertv et al.l (pOM ). 

We ignore direct interactions of CRs with the ambient 
cluster gas, e.g., Coulomb interactions, hadronic colli- 
sions and interactions through the generation of hydro- 
magnetic waves, which may become important especially 
when the cavities are disrupted and wh en the CRs are 
mixed with the ICM (|Guo fc Ohll2008l ). These inter- 
actions may provide significant heating effects for the 
thermal gas. One byproduct of these interactions is the 
generation of 7-rays through the decay of neutral pions, 
which ma y be studied by the recently- launche d Fermi 
telescope (jAndo fc Nagaf 120081: IMathewsj 12009ft . Previ- 
ous cluster observations with the High Energy Stereo- 
scopic System (HESS) have not detected these 7-ray sig- 
nals in the VHE (> 100 GeV) band, providing upper 
limits on 7-ray fluxes and ratios of cosm ic ray pressure to 
therm al pressure in galaxy clusters (e.g.. lDomainko et al.l 
l2009h . 

The center of each cavity in the cluster MS0735 is lo- 
cated at a projected cluster-centric radius of around 170 
kpc, which is much larger than the average projected 
cluster radius of 2 kpc for X-ray cav ities typically ob- 
served in clusters ( Bir zan et al.l 120041 ) . If strong AGN 
outbursts having energies similar to that of MS0735 ere- 
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ate X-ray cavities much closer to the cluster center than 
those observed in MS0735, they produce much stronger 
heating in the cluster cool core and may even fully de- 
stroy it, tra nsforming a cool core cluster to a non-cool 
core cluster (|Guo &: Ohll2009( ). It is of great interest to 
investigate if metallicity peaks at cluster centers are also 
destroyed in this process, producing much flatter metal- 
licity profiles as observed in non-cool core syst ems (e.g., 
IDe Grandi fc Molendil l200l IBaldi et alJl200l . We will 



thoroughly investigate these models in a follow-up paper. 



We thank Fabrizio Brighenti for useful discussions, and 
the anonymous referee for a positive and helpful report. 
Studies of the evolution of hot cluster gas at UC Santa 
Cruz are supported by NSF and NASA grants for which 
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